Topological Vortex Formation in BEC under Gravitational Field 
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Topological phase imprinting is a unique technique for vortex formation in a Bose- Einstein conden- 
sate (BEC) of alkali metal gas, in that it does not involve rotation: BEC is trapped in a quadrupole 
field with a uniform bias field which is reversed adiabatically leading to vortex formation at the 
center of the magnetic trap. The scenario has been experimentally verified by MIT group employ- 
ing 23 Na atoms. Recently similar experiments have been conducted at Kyoto University, in which 
BEC of 87 Rb atoms has been used. In the latter experiments they found that the fine-tuning of 
the field reverse time T rev is required to achieve stable vortex formation. Otherwise, they often 
observed vortex fragmentations or a condensate without a vortex. It is shown in this paper that this 
behavior is attributed to the heavy mass of the Rb atom. The confining potential, which depends 
on the eigenvalue ms of the hyperfine spin F along the magnetic field, is now shifted by the grav- 
itational field perpendicular to the vortex line. Then the positions of two weak-field-seeking states 
with tub = 1 and 2 deviate from each other. This effect is more prominent for BEC with a heavy 
atomic mass, for which the deviation is greater and, moreover, the Thomas-Fermi radius is smaller. 
We found, by solving the Gross-Pitaevskii equation numerically, that two condensates interact in a 
very complicated way leading to fragmentation of vortices, unless T rcv is properly tuned. 

PACS numbers: 03.75.Mn, 03.75.Kk, 67.57.Fg 
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I. INTRODUCTION 

Alkali atom gas undergoes phase transition to a super- 
fluid phase at very low temperatures (~ 0.1 mK) to form 
a Bose- Einstein condensate (BEC) [H 0, 0, ■ This new 
BEC is different from the conventional superfluid 4 He 
in that (i) it is a weakly coupled system for which the 
Gross-Pitaevskii equation is applicable with high preci- 
sion at low temperatures (ii) alkali atoms have hyperfine 
spin degrees of freedom F and, accordingly, the order 
parameter has 2F + 1 components 0, while 4 He atoms 
have no such internal structure (iii) external magnetic 
field couples with the hyperfine spin, allowing for easy 
control of the condensate order parameter, among oth- 
ers. For example, the atoms 23 Na and 87 Rb have internal 
states F = 1 and F = 2, for which the number of com- 
ponents is 3 and 5, respectively. 

Associated with BEC is a quantized vortex, which 
demonstrates the existence of coherent property and su- 
perfluidity of the condensate. The study of vortices in 
BEC of alkali atoms started as soon as the BEC was 
discovered and several scenarios for vortex formation 
were proposed, namely (i) dynamical phase imprinting 
El i (ii) o ptic al spoon 0, O El (iii) rotating asym- 
metric t rap 113,113 and (iv) topological phase imprinting 
LUl LLj- UH flTl Hq . We will be concerned in the present 
paper with a vortex formation in the last scenario under 
the gravitational field. 

In topological phase imprinting, the BEC is confined 
in a quadrupole magnetic field and a bias field orthogonal 
to it. As the bias field is adiabatically reversed, the con- 
densate spin follows the direction of the local magnetic 



field and position dependent order parameter phase is 
imprinted in BEC in such a way that a vortex is created 
along the axis of the trap when the bias field is reversed. 
This scenario has been verified in beautiful experiments 
at MIT and reported in Leanhard et. al. [18], in which 
F = 1 and F — 2 hyperfine states of 23 Na have been 
utilized. The analysis based on Berr y's phase llfjlljjj as 
well as numerical solutions 0, 0, , predicts that a 
vortex thus created has the winding number IF ', which 
was also verified in Ha. 



In recent experiments conducted at Kyoto University 
[2fj| . where BEC of 87 Rb has been used, they observed 
that fine-tuning of the field reverse time T rcv is required 
for topological formation of a stable vortex. Vortex frag- 
mentation or no vortex has been observed for different 
choices of T rov . We show, in this paper, the existence 
of a narrow window in T rov for successful formation of 
a stable vortex is attributed to a subtle role played by 
the gravitational field and the heavier mass of Rb com- 
pared to 23 Na. 



This paper is organized as follows. In Sec. II, we out- 
line the formalism of the problem based on the multi- 
component spinor Gross-Pitaevskii equation. Section III 
is devoted to the numerical results which verify our state- 
ments. It is shown that our T rcv for stable vortex forma- 
tion approximately reproduces that obtained in Kyoto 
experiments. Conclusions and discussion are given in Sec. 
IV. 
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II. TOPOLOGICAL PHASE IMPRINTING 
UNDER GRAVITATIONAL FIELD 

We will be mostly concerned with BEC of 87 Rb in F = 
2 hyperfine state. The case F = 1 may be analyzed in 
a similar manner with considerably less effort. We first 
introduce the order parameter for F = 2 BEC and the 
Gross-Pitaevskii equation which describes the dynamics 
of the condensate |2lL |22| . Topological vortex formation 
in the presence of a gravitational field is then discussed 
in detail. 



where the last equality defines \2)b- This state is under 
an attractive force toward a region with smaller B = \B\ 
and hence is called a weak-field-seeking state (WFSS). 
We also obtain \1)b starting from |1) in Eq. QJ, which is 
also called a weak-field-seeking state. The states | — 1)b 
and | — 2}b, obtained from | — 1) and | — 2), respec- 
tively, are called strong-field-seeking states (SFSS) while 
the state \0)b is independent of the magnetic field and 
will be called the neutral state (NS). 



B. Gross-Pitaevskii Equation 



A. Order Parameter 

Let us consider a uniform BEC with hyperspin F = 2. 
The internal state of the condensate may be classified 
according to the eigenvalue m z of the operator F z , where 
— 2 < m z < 2. The corresponding eigenvector is denoted 
as \m z ): F z \m z ) = m z \m z ). As for the eigenvectors we 
employ the following convention 

|2> = (1,0,0,0,0) T , |1) = (0,1,0,0,0) T , 

|0> = (0, 0, 1, 0, 0) T , I - 1) = (0, 0, 0, 1, Of, (1) 

|-2) = (0,0,0,0,1) T , 

where T stands for transpose operation. An arbitrary 
order parameter is then expanded as 



I*) 



2 

E 

rn z — — 2 



= (* 2 ,*i,*o,*-i,*- 2 ) T . (2) 



In case the external magnetic field B is position- 
dependent, the classification of the hyperfine state in 
terms of m z is not appropriate and, instead, the eigen- 
states of Fb = F ■ B / B must be utilized. 

Consider the order parameter 



|*)o = /o|2), 



(3) 



which satisfies the eigenvalue equation F z \^)q = 2|\E f ) - 
The corresponding eigenstate with respect to an arbi- 
trary F B , F B \W) = 21*), with 



sin /? cos a 
B(r) = B \ sin (3 sin a 

cos (3 



(4) 



The dynamics of F = 2 BEC in the limit of zero tem- 
perature is given, within the mean field approximation, 
by the time-dependent Gross-Pitaevskii (GP) equation 
which is extended so that the spin degrees of freedom are 
taken into account as |2l| 
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where summations are over k = x,y, z and —2<n,p,q< 
2. Here, M is the mass of the atom and the gravita- 
tional force is assumed to be along the positive x axis. 
The Larmor frequency is defined as hu>L k — where 
7 M ~ /is is the gyromagnetic ratio of the atom, /is being 
the Bohr magneton. The interaction parameters {gi} are 
expressed in terms of the s-wave scattering length clf, F 
being the total hyperfine spin of the two-body scattering 
state, and are given by 
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The parameters are oq = 4.73 nm, 02 
04 = 5.61 nm for 87 Rb atoms. 



5.00 nm and 



is obtained by rotating l^g) through Euler angles a, (3 
and 7. Explicit calculation shows that 

I*) = exp(-mF z )exp(-z/3F y )exp(-i7F z )|* ) 



= foe 



-217 



COS 4 § \ 



2e- la cos 3 f sin f 
V6 cos 2| sin 2| 

2e ia cos f sin 3 f 
e 2la sin 4 f 



= h\2) B , (5) 



C. Topological Phase Imprinting 

Suppose a WFSS with raj = 2 is placed in a Ioffc- 
Pritchard trap. We consider a condensate uniform along 
the z axis and confined within the xy plane by the 
quadrupole magnetic field 



B±(r) 



B±(r) cos(- 
B ± (r) sin(- 
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where B±(r) is, to a good approximation, linear in r 
within the Thomas- Fermi radius; B±(r) = B'r. Here 
(r, (j>, z) are the cylindrical coordinates where the axis of 
the magnetic trap is taken as the z axis. To prevent the 
condensate from escaping from the trap through Majo- 
rana flips, a uniform bias field B z (t) must be applied 
along the z axis, where we have explicitly written the t 
dependence of the bias field, which is necessary for topo- 
logical phase imprinting. The total magnetic field is thus 
given by 



B(r,t) = B ± (r) + B z (t) = 



B±(r) cos(- 
B±(r) sin(- 
B z (t) 



(8) 



The condensate would be cylindrically symmetric around 
the z axis if it were not for the gravitation field. In the 
presence of the gravitational field with a potential —Mgx, 
in contrast, the condensate is not symmetric any more. 

The initial condensate profile is determined by solving 
the stationary GP equation, in which the ms = 2 WFSS 
|J5J is substituted. The angles in JSJ) are now given by 



a 



[3 = tan" 



= — tan 
-i B ± (r] 



-i y 

X 



B z (t) 
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(9) 



It is assumed that the initial state is vortex free, which 
forces us to choose a — —'){= — 0) in Eq. (J5J. The 
ground state condensate profile corresponds to the lowest 
energy eigenvalue of the equation 



2M 



where v r , 



-Mgx (f v m ) + (gi + 4g 2 )f^v m 

+fujj L f v m = fJ>fov m , (10) 
= (m|2) B , V 2 = d 2 /dx 2 + d' 2 /dy 2 and 

hiu L = 1 ^B' 2 {x 2 +y 2 ) + B z {t) 2 . (11) 



The eigenvalue [i is identified with the chemical potential. 
If {w m } t is multiplied from the left and J2 m \ v m\ 2 = 1 and 
similar identities are employed, we obtain the reduced GP 
equation for fa; 



1L 

2M 



V 2 / + J2 v *n {(V 2 « m )/o + 2Vv m ■ V/ } 
Mgxf + ( 9l + 4 ff2 )/ 3 + tkJ L f = nf , (12) 



The initial condensate profile and the chemical potential 
H are obtained by solving (|T21) numerically. 
The uniform bias field is then reversed as 



BS) 



BM 1 



2t 



T 



-B z (0) 



0<t< T rcv 



(13) 



It has been shown, in the previous works |16l 1 1 7| with- 
out taking gravity into account, that a vortex formation 
takes place at the center of the condensate when the 
bias field has been reversed. The bias field disappears 
at t ~ T rcv /2 and the tub — 2 condensate transforms 
into —2 < nifl < 1 components through Majorana flips. 
The components with ms < are not confined and es- 
cape from the trap. The final condensate is then made 
of a mixture of the ms = 1 and = 2 WFSS. It is 
found from Eq. (JSJ) with (3 = n that the vortex in the 
m.B = 2 component has the winding number four while 
a similar condensate amplitude for = 1 tells us that 
the vortex in this component has the winding number 
three ^ is essential, therefore, to employ the full 

spinor GP equation @ to analyze this system. 

In the presence of the gravitational field, however, the 
center of the condensate shifts from the axis of the mag- 
netic trap and, moreover, the minimum of the trapping 
potential, including gravity, depends on the hypcrfinc 
state due to the last term in © . To be more concrete, the 
ms component of the condensate in the Ioffe-Pritchard 
trap is subject to the external potential 

V(x,y,t) = ^m B j^B' 2 (x 2 + y 2 ) + B z (t) 2 - Mgx. 

(14) 

Note that the center of the quadrupole field is (0, 0) while 
the minimum of the potential is 



(x ,y ) = 



2Mg\B z (t)\ 



B'JB' 2 m 2 B j 2 -4M 2 Q 2 



(15) 



where we have assumed that B'tubJ^, > 2Mg, which 
always holds in actual experiments. Otherwise, there is 
no stable minima in the x-direction and the condensate 
cannot be confined. Observe that the deviation xq has 
the same time dependence as |£?^(t) | and that it is larger 
for tub — 1 than for raj = 2, implying that there is a 
finite distance between the centers of the tub — 1 and 
tub = 2 condensates. It is important to note also that 
xq is proportional to the atomic mass M, provided that 
B'm_B7 M 3* 2-Mg, and hence the effect of the gravity on 
the condensate is more prominent for atoms with heavy 
mass. 

The potential V (x, y, t) is approximated by a displaced 
harmonic potential when \B Z \ 3> |-B_l|, taking the form 



V(x,y,t) 



Z\B z (t)\ + y2) ~ MgX + 2 mBl ^ B ' M 



Mg 2 



"HO 



(*) 



2o4 (i) ' 2 



-m S 7 M |S z (t)|, 



where 



w HoW 



2M\B z (t)l 



(16) 



(17) 
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The displacement xq = g/oj^ (t) found from l|16|) is in 
agreement with the result (|15fl when B'itibj^ ^S> 2Mg. 
The frequency who = w ho (0) serves as an energy scale of 
the system. The gravitational potential in Eq. (|14|l leads 
to two remarkable effects on vortex formation. First, for 
F = 2 BEC, there are two WFSS, one (m B = 1) has 
three units of vorticities while the other (tub = 2) has 
four units. Their positions deviate due to the difference 
in rris, which makes the interaction between the con- 
densate even more complicated. Secondly, the minimum 
of the potential shifts from the center of the quadrupole 
field, where the topological vortex formation takes place. 
A highly quantized vortex with winding number \n\ > 2 
is unstable against decay into \n\ singly quantized vor- 
tices. The lifetime of a highly quantized vortex is longer 
if its axis is at the center of the condensate. If it is 
off-centered, in contrast, it is expected that the vortex 
moves toward the edge of the condensate and split into 
\n\ singly quantized vortices. Therefore a topologically 
created vortex under gravity will be unstable against de- 
cay unless some caution is exercised, see the next section. 
If the deviation of the center of the quadrupole field is 
large compared to the Thomas-Fermi radius, the "vortex 
formation" takes place outside the condensate and it is 
impossible to observe the vortex. 

These effects are more prominent in 87 Rb compared 
to 23 Na. In Kyoto experiments, the parameter who — 
2tt x 330 [rad/s] is taken to be the same as the MIT 
experiment. Then the Thomas-Fermi radius ttf = 
fl HO\/2/VA/w^ has mass dependence ~ M -1 / 2 through 
ano = \J fi/MuHO since the factor /i/M roughly takes 
the same value for both experiments. 

The above consideration shows that the behavior of the 
condensate depends heavily on the atomic mass M and 
the reverse time T rev - We solved the spinor GP equation 
© numerically for various T rcv for both 23 Na and 87 Rb. 
In the next section, we outline mainly the results for the 
F = 2 condensate of 87 Rb and then compare them with 
those for 23 Na. 



III. NUMERICAL RESULTS 

We have solved the stationary GP equation l|12(l nu- 
merically to obtain the vortex-free initial condensate pro- 
file fo(x,y) and the chemical potential /i. We then 
solved the time-dependent spinor GP equation JSJ) with 
^ m = fov m as the initial state. We had to pay special 
attention to particles escaping form the trap. This is 
taken care of by introducing an appropriate absorption 
potential beyond which the condensate amplitude decays 
exponentially. We found slight interference patters due 
to reflection at the absorption potential, although we be- 
lieve that this does not cause any qualitative change in 
conclusion. 

As for numerical values, we take the following numbers 
reproducing the Kyoto experiments: N = 5x 10 5 , who — 
2tt x 330 [Hz], uj z = 2tt x 17 [Hz]. The Thomas-Fermi 



radius is then ttf = 3.1 /im. The initial bias field is 
B z (0) = 0.5 G, which is then reduced to B Z (T ICV ) — 
-0.5 G. 

It turns out that a stable vortex formation takes place 
only for T rev in a small window around T rcv ~ 2 ms. The 
behavior of the BEC is qualitatively different for T rcv 
out of this window from the case T rcv ~ 2 ms. This is in 
harmony with the observation made by the Kyoto group 
1201 . Therefore 4 ms, 2 ms and 1.5 ms as 

typical values representing the three classes of T rev . 

A. T rov = 4 ms 

The change in B z is slow in this case and each WFSS 
component of the condensate follows its equilibrium po- 
sition l|15f) as B z is reversed. There is essentially a single 
condensate with ms = 2 at t < T rcv /2. Majorana flips 
take place at t ~ T rcv /2 in the vicinity of the point where 
the magnetic field vanishes. The center of tub = 2 con- 
densate is separated from that of ms = 1 condensate 
from then on. There exists a vortex with the winding 
number four (three) in the raj = 2 (ms = 1) conden- 
sate, which are created at the center of the quadrapole 
field at around t ~ T rev /2. The center of the condensate 
is at the center of the quadrupole field only at t = T Tev /2 
and its subsequent position follows the minimum of the 
potential, which depends on tub- Then these compo- 
nents interact in a complicated manner with each other 
and vortex fragmentation takes place. 

Figure 1 shows the total density profile J2 m l^ml 2 and 
the component profile \^m B 1 2 f° r m B — 1,2 with T lcv = 
4 ms. 



B. T rov = 2 ms 

A window for stable vortex formation exists in the 
vicinity of T rev ~ 2 ms. The condensate is pushed up by 
reducing B z (t) and the center of the condensate reaches 
the center of the quadrupole field at t ~ T rev /2 as in 
the previous case. The minimum of the potential fol- 
lows the path it took for < t < T rcv /2 backward af- 
ter t = Trcv/2. The condensate, however, fails to follow 
the potential minimum since it has been accelerated up- 
ward. During the acceleration, all atoms are in a single 
hyperfine state and equally accelerated. Following the 
Majorana flips at t ~ T rcv /2, each component eventu- 
ally reaches a peak with some delay, after which it starts 
to fall chasing for the corresponding potential minimum. 
The delay is longer if the acceleration is larger, i.e. T rcv 
is smaller, and the strength of the confinement is weaker. 
The latter condition depends on tub as shown explicitly 
in Eq. fTfl . 

If is fine-tuned to be around 2 ms, the distance 
between the two condensates with = 1 and ms = 2 
happens to be small at t = T rev in spite of separated 
potential minima for respective components. Therefore 
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the iriB = 2 vortex with the winding number four is not 
yet disturbed by the ijib = 1 vortex with the winding 
number three, see Fig. 2. Observe that the center of the 
itib = 1 component is below (above) that of the ms = 2 
component at t = T rcv in Fig. 1 (Fig. 3). Vortex frag- 
mentation occurs at t > T rcv when the distance between 
the condensates becomes larger as in the previous case. 

C. T rcv = 1.5 ms 

Since the change in B z (t) is fast, the condensate can- 
not follow the minimum of the potential (|15f) and does 
not reach the center of the quadrupole field at t = T rcv /2. 
The center of the vortex always appears in the center of 
the quadrupole field and thus the vortex axis appears 
above the center of the condensate in this case. The 
vortex thus created has the winding number four and a 
highly quantized off-centered vortex is expected to decay 
quickly into four singly quantized vortices |23|. This is 
in fact manifest in Fig. 3 at t — T rov . In Kyoto experi- 
ments |2JJ, the wavy edge of the condensate is observed 
when T rov is smaller compared to that for stable vortex 
formation. This might be attributed to the fragmented 
vortices. 



IV. CONCLUSIONS AND DISCUSSION 

In summary, we have shown that the effect of the grav- 
itational field on the topological vortex formation in BEC 
of heavy atoms, such as 87 Rb, is not negligible. Unless 
the reverse time T lcv is fine-tuned, fragmentation of a 
multiply quantized vortex takes place during topological 
phase engineering. This fragmentation is attributed to 
the existence of two WFSSs with tub — 1 and 2, whose 
centers are separated, under gravity, by distance compa- 
rable to the Thomas-Fermi radius. WFSS with ms = 1 
supports a vortex with three units of winding number 
while that with tub = 2 supports four units. 

When the two vortices have the same center, there 
exists an axially symmetric stable solution ^tJ- If t ne 
vortex centers are displaced from each other, however, 



the two WFSSs interfere with each other through the g?- 
and (?3-terms in Eq. © so that the deviation from axial 
symmetry is further amplified, leading to fragmentation 
observed in [2(j ■ 

When the field is reversed so slowly that each hypcrfine 
state follows its potential minimum, the center of two 
hyperfine components deviate from each other as soon 
as Majorana flips take place at t ~ T rcv /2, leading to 
vortex fragmentation. When the field is reversed so fast, 
in contrast, the condensate fails to follow the potential 
minimum and the deviation of two components at t = 
T rev is minute. The vortex is, however, imprinted near 
the edge of the condensate, which leads to spliting of the 
multiply quantized vortex into singly quantized vortices. 
The fragmentation pattern often has two-fold symmetry, 
which originates from the vertical displacement of two 
condensates |23l | . To observe a stable multiply quantized 
vortex, the field reverse time must be taken somewhere 
between above two cases. 

We have also studied the effect of the gravitational field 
on the BEC of 23 Na with F = 2, which corresponds to 
the MIT experiments. We have shown that there exist 
a similar window of T rGV for a topological formation of a 
stable vortex in this case. 
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FIG. 1: (Color) Time development of the condensate with T ICV — 4 ms. The white lines show the coordinate system, the 
origin of which is the center of the quadrupole field where the vortex formation takes place. Upper panel shows that total 
density of the condensate. The number in the panel shows the normalized time t/T Tev . The middle (lower) panel the density 
of ms = 2 (ms = 1) component. The cross (x) shows the minimum of the potential for the respective component. Observe 
that the center of each ms component follows the minimum of the respective potential and accordingly the center of the the 
ms = 1 component is below that of the tub = 2 component at t = T rcv . 
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FIG. 2: (Color) Time development of the condensate with T lev = 2 ms. See the caption of Fig. 1 for the definition of the 
symbols. The condensate is pushed up while the potential is deformed by reversing B z . The ms = 1 component is still above 
the minimum of its potential at t — T rcv so that the center of the component is close to that of the ms = 2 component. 




FIG. 3: (Color) Time development of the condensate with T rcv = 1.5 ms. See the caption of Fig. 1 for the definition of the 
symbols. The condensate is pushed up by a quick deformation of the potential and the centers of both ms components are still 
above the minima of the potentials at t — T ICV . This is more eminent for ms = 1 component so that the center of the ms = 1 
component is above that of the ms = 2 component. 



